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Matrix elements of heavy-light flavor currents play an important role in modern particle physics 
and precise theory predictions are of interest for phenomenology. Heavy Quark Effective Theory 
(HQET) is a valuable tool to obtain such predictions. In the HQET matching program of the AL- 
PHA collaboration presently only the temporal component of the axial vector current is included. 
Extending the matching to the temporal component of the vector current and the spatial compo- 
nents of the axial vector current thus seems desirable. Here we present a recent one-loop study in 
lattice perturbation theory to test two candidate matching observables for these currents for their 
quality to guide future non-perturbative investigations [1]. 
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1. Introduction 

Matrix elements of heavy-light flavor currents are important theoretical input to constrain 
physics of and beyond the standard model through decays of heavy mesons. A theoretically clean 
way to treat systems with one single heavy quark is Heavy Quark Effective Theory (see [2] and 
references therein), short HQET. From the theory point of view, especially leptonic decays of heavy 
mesons are relatively easy to treat. Hence next to the parameters of the action, the HQET program 
of the ALPHA collaboration [3-8] included the temporal component of the heavy-light axial vector 
current form the very beginning. 

Because of the current big interest in flavor physics within the community, it is a logical next 
step to include all components of the vector and axial vector current to make theory predictions 
for matrix elements such as for example the one entering in the exclusive semi-leptonic decay 
B -)• nlv. 

1.1 HQET and matching 

HQET is a systematic expansion of QCD with one heavy quark in inverse powers of its mass. 
In addition to the parameters in the HQET action, one has to take care of the coefficients of the 
currents under investigation as well. The inclusion of all components of the vector and axial vector 
current increases the numbers of HQET parameters to be fixed at order 1/mtoa total of 19. Con- 
sidering this rather big number of matching conditions, it seems useful to test them first in relatively 
cheap perturbative calculations. This possibility is even more attractive since we have the pastor 
software package [9] for automated lattice perturbation theory calculations at our disposal. Even 
though perturbative one-loop calculations performed with pastor only require a moderate effort 
we do not intend to perform the matching of all 19 parameter here. Instead we will focus on the 
renormalization constants of the temporal component of the vector and the spatial components of 
the axial vector current. These play an important role since as stated above they were not included 
in the matching by the ALPHA collaboration so far and enter already at the lowest order in HQET. 

The matching procedure consists basically of evaluating the same number of observables <J>, 
as the number of parameters to be fixed - in QCD and HQET - and setting <J>^ QET = <1>P CD to 
determine (in fact define) the HQET parameters. The main objective of this work will be to in- 
vestigate possible matching conditions for the currents mentioned above and estimate the size of 
higher order contributions in the 1/m-expansion. For a precise determination of the parameters, 
one certainly wants these higher order effects to play as little a role as possible. 

2. Matching of flavor currents 

We work in the Schrodinger functional [10, 1 1] with two mass-less quark flavors, m\=m2= 0, 
and a third heavy one whose renormalized mass is as usual given by the dimension-less parameter 
z = m^L = mL. The boundary fields are given by Q, and the phase angle for all three flavors is 
equal and denoted 6. We are interested in the temporal component of the vector current and the 
spatial components of the axial vector current of a heavy and a light quark, 

M x ) = W3(x)YkY5V2(x), V Q (x) = \j/ 3 (x)Yoy 2 (x), (2.1) 
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or to be more precise their matrix elements between heavy-light pseudoscalar states. To this end, 
we define in a first step the three-point functions (with implicit improvements b A ,by [12]) 

Z X L X L X 

f?°(xo;e,z) = -^E(Ci)sC3AoWC2)sCi'>- (2-2) 

Z X 

To define observables that are suitable to extract the desired matrix elements [13], we form the 
ratios 

V [/i(M/i(e,o)] 1/2 y V IMM^M)] 172 / 

/i(e,z) = -f KCirsCsCsrsCO, *i(M =-f KCinCaCsrid'}- (2.4) 

X X 

The renormalization constants Zv,Za are fixed by Ward identities [14, 15]. The quark masses are 
renormalized using the lattice minimal subtraction scheme [16, 17] at one-loop order, 

m = Z mM (gl,alU,)m q [l+ab m (gl)m q ], m q =m -m c , (2.5) 

b m = -0.5- 0.07217(2) C f g%, Z mSat (gl,a^) = 1 - _l og (ajt)gg. (2.6) 

27T Z 

2.1 Static approximation 

We can use the symmetries of the static theory to establish a relation between the renormalized, 
renormalization group invariant (RGI) axial vector and vector currents [2], 

y^HQET = Cps (M b /A m )Z^ RGl (g )^ A (go)Vr, (2-7) 
A?* ET = Cy(M b /A m )Zri G MZ^ A (go)Ar- (2-8) 
At one loop level, we have in the lattice minimal subtraction scheme 

(C)o(m) =zZMz^ A vt\ <AT)M=z&MZftAr, (2-9) 

= 1 " X>log(aM)s§ + 0(*g), 4% = 1 + (^a) (1 ^o + 0(^), (2.10) 
with (c.f. [18-20]) 

Z^ A ) (1) =0.0521(1), 7o = -^- (2.H) 
We define the static counterparts to <1> A * and <J> V °, 

x Ak (n) = log (z&ooz^) +< re , x Vo (m) = log (zSStOijz^k) +*vT> ( 2 - 12 ) 

with and Xy^ re given (without the renormalization factors) as in (2.3) with the heavy quark 
flavor replaced by a static one. We may then expect that at one loop level 

O v »(L,m) =Bf t gl+X Y (m) + 0(l/m), ^(L,m)=Bf x gl+X A (m) + 0(\/m), (2.13) 

where we have from [2], using the ratio of Cy and Cps, 

Bf l = -0.137(1), Bf l -Bf l = 0.016900. (2.14) 
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3. Results 

3.1 Tree level 

Tree level results for <I> Vo and <I> A * are shown in figure 1, the points at 1/z = correspond 
to the static values and Xy Q . For 6 = 0.5, 1.0, the dependence on 1/z is linear with a small 
slope. In the interesting region of z = 10, which is the typical matching point for B physics used 
in [21], the 1/z corrections at tree level are of a few percent. For 6=0, there is no dependence on 
z and L/a at tree level for both observables. The continuum limit at tree level was extracted from 
computations with L/a up to 200, using the fitting procedure explained in [22]. 

3.2 One loop 

The observables fx, <J> V °, and <J> A * were calculated at one loop level for z = 4,6,8, 10 (and at 
z = for f\, kx), with lattice resolutions up to L/a = 40, and for 6 G {0.0,0.5, 1.0}. Furthermore, 
we evaluated the static counterparts ff m , <E> v o> stat and A * stat for lattices with L/a up to 28. No 
bigger lattice sizes are required for the HQET quantities, since their continuum limit is easier to 
obtain due to a weaker a/L-dependence. The continuum values are presented in figure 2, the points 
for 1/z = are again the static values on which we will comment on shortly. The extraction of 
the continuum values is done with the method mentioned above. To establish the connection to the 
static limit, we define the one loop quantities 

GL 1) (z)=^ 1 )(z) + 7olog(z)^>4 1) , ae{V ,Ax}, (3-1) 
< =4r (1) + (^%) (1) +sr-7olog(«/L), (3.2) 
H (i) =x b.e, ( i) + ^« +Br _ 7olog(a/L) (33) 

The fit functions included in the plots are of the form G« (z) = + /z + c" /z 2 - Even though 
they are only indicative (since possible terms of the form log (z)/z" are not included), one can 
anticipate that the higher order collections corresponding to C2 are rather small. We did not attempt 
to fit the logarithmic terms due to the small number of available data points. 

In the case of <J> A *, a value of 8 = seems to minimize the higher order corrections in 1/z both 
at tree level and at one loop. For <J> V ° , the value 0=0.5 seems to be a good compromise considering 
the corrections at the two orders in perturbation theory we investigated. One should always keep 
in mind that the small observed dependence on 1/z will be eliminated once the 1/wih corrections 
are included in the effective theory. Finally only the l/z 2 -terms, manifesting themselves in the 
curvature of the fits, will remain as corrections. 

4. Conclusions 

We have seen that the proposed observables for the matching are dominated by the lowest 
orders in l/m and hence one would strongly suggest to use them in a non-perturbative study. For 
the phase angle 6, we would suggest the value 0.5 as a compromise to minimize 1/z-effects at 
one-loop and tree level, but the current study does not indicate a strong preference for this value. 
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A further investigation including the matching conditions for all 19 parameters is currently carried 
out within the ALPHA collaboration. 

We could also demonstrate the usefulness of the pastor software package for applications 
beyond the classical realms of lattice perturbation theory such as renormalization and improvement. 
The publication of an in-depth description of pastor along with its source code is planned for the 
near future. 
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A. Plots 
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Figure 1: <I> Aa and <J> V ° at tree level in the continuum limit. The point size is bigger than the error. 
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Figure 2: and Gyj in the continuum limit. The dashed lines are a linear fit in l/z, setting C2 = 0. The 
data points at z = 4 are not included in the fits. 
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